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We show that a large class of two-field models of single-bubble open inflation do not lead to infinite 
open universes, as it was previously thought, but to an ensemble of very large but finite inflating 
'islands'. The reason is that the quantum tunneling responsible for the nucleation of the bubble 
does not occur simultaneously along both field directions and equal-time hypersurfaces in the open 
universe are not synchronized with equal-density or fixed-field hypersurfaces. The most probable 
tunneling trajectory corresponds to a zero value of the inflaton fleld; large values, necessary for the 
second period of inflation inside the bubble, only arise as localized fluctuations. The interior of each 
nucleated bubble will contain an infinite number of such inflating regions of comoving size of order 
^ ■ 7~^) where 7 depends on the parameters of the model. Each one of these islands will be a quasi-open 

O ' universe. Since the volume of the hyperboloid is infinite, infiating islands with all possible values 

, of the field at their center will be realized inside of a single bubble. We may happen to live in one 

of those patches of comoving size d < 7~^, where the universe appears to be open. In particular, 
we consider the "supernatural" model proposed by Linde and Mezhlumian. There, an approximate 
C/(l) symmetry is broken by a tunneling field in a first order phase transition, and slow-roll inflation 
^ ' inside the nucleated bubble is driven by the pseudo-Goldstone field. We find that the excitations 

' of the pseudo-Goldstone produced by the nucleation and subsequent expansion of the bubble place 

severe constraints on this model. We also discuss the coupled and uncoupled two-field models. 

(N 
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I. INTRODUCTION 



^ . In models of open inflation, which lead to a density parameter f^o < 1, the "horizon" and "flatness" problems 
' are solved by two very different mechanisms. Ahhough open inflation can be reahzed with a single scalar, reahstic 
, models look more natural when the task of solving each one of these two problems is entrusted to a different scalar 
field. Nevertheless, models with two fields introduce a host of new effects which should be carefully investigated. In 
particular, as we shall see in this paper, most of the two-field models that have been recently proposed do not give 
' rise to an infinite open universe but to a large inflating island of finite size: a quasi-open universe. 
\ The picture of open inflation is the following. The universe starts in a de Sitter phase driven by the potential 
energy of a scalar fleld a which is trapped in a false vaccuum. This false vacuum decays through quantum tunneling, 
and spherical bubbles of true vacuum nucleate in the smooth de Sitter background. After nucleation, the bubbles 
expand with constant acceleration, following a "trajectory" which is invariant under Lorentz transformations 0(3, 1), 
see Refs. and Fig. [l]. Since this is also the symmetry of an open Friedmann-Robertson- Walker (FRW) universe, 
the <T — const, surfaces in the interior of the bubble can be identified with the t = const, sections of an open universe 
In this way the symmetry of the bubble takes care of the "homogeneity" problem. A second period of 'slow-roll' 
inflation inside the bubble, lasting for approximately 60 e-foldings, would solve the "flatness" problem [§-0. 

In open (and in quasi-open) inflation, the dynamics of bubble nucleation and subsequent expansion turns out to be 
very important in determining the spectrum of gravity waves and density perturbations. The reason is that, unlike 
the case of stardard inflation, the amount of slow-roll inflation is minimal and the "initial conditions" right after 
the bubble nucleates are not washed out completely. Thus, for instance, quantum fluctuations of the slow-roll field 
generated outside the bubble can penetrate to the interior g] , causing perturbations whose wavelength is larger than 
the curvature scale. These are the so-called supercurvature modes Also, the scattering of tensor modes off the 
bubble wall determines the spectrum of very long wavelength gravitational waves JTo| , pd] |. In the limit of a weakly 
gravitating wall, this effect can be alternatively described as a fluctuation of the bubble wall itself, which induces 
supercurvature anisotropics inside the bubble JT^-^lSt. 

In principle, tunneling and slow-roll can be done by the same scalar field Q|, but this requires a very special form 
of the inflaton potential V . Denoting hy H = (SttGV^/S)^/^ the Hubble rate during inflation, a sharp barrier where 
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^ p = const- 



FIG. 1. Conformal diagram showing a bubble expanding in a de Sitter background. The bubble wall is represented as 
the thick grey line starting at the point P. It expands with constant acceleration along a. rj — const, surface. The FRW open 
universe is inside the future light cone from the point O, which is the center of symmetry of the bubble solution. 

V" ^ is necessary^ for bubble nucleation But this barrier must be right next to a flatter region where 

V" ^ H^, which is needed to make slow-roll possible. Moreover, the duration of slow- roll, which depends on the 
length of the plateau in the inflaton potential, has to be fine-tuned to some extent, because a few e-foldings more or 
less can make the difference between an almost flat universe and an almost empty one. 

As mentioned above, models with two fields were introduced in order to overcome these difficulties; one doing the 
tunneling and the other doing the slow-roll fn this way, the coexistence of two different mass scales seems more 
natural. Also, it was argued that in some models the value of the slow-roll field after bubble nucleation can be different 
in each nucleated bubble, and hence the duration of open inflation would be different in each one. As a result, for 
a given temperature of the CMB, one would obtain a different value of the density parameter in each universe, and 
there would always be some open universes with a density parameter in the interesting range . 

The purpose of this paper is to show that in models of this sort, with variable fi, the picture is actually more 
complicated. Indeed, instead of an infinite open universe inside of each bubble, what we find is an infinite number of 
infiating islands of finite size inside each bubble. 

Quasi-open universes are not entirely new. The simplest two-field model of open infiation, where the tunneling 
field a and the slow- roll field (p are decoupled, is actually a quasi-open one, as emphasized by the authors of Rcf . \p \. 
Quasi-openness is in principle not a desirable feature, since to a typical observer, the universe looks anisotropic [[l9 |. 
In the simple "decoupled" model, this "classical" anisotropy is large and, combined with the effect of quantum 
"supercurvature" fiuctuations mentioned above it basically rules out the model [ p^ . 

To circumvent this problem, Linde and Mezhlumian introduced a class of two-field models where the slow-roll field 
is coupled to the tunneling field. As we shall see, these models are also quasi-open. This does not mean that they 
are not good cosmological models. If the co-moving size of the inflating islands is sufficiently large, then the resulting 
classical anisotropy may be unobservable. Even so, the fact that these islands are finite leads to a dramatically 
different picture of the large scale structure of the universe in open models. 

The plan of the paper is the following. In Section II we briefly review the "supernatural" model of inflation, 
introduced in Ref. ||^. In this model, an approximate U{1) symmetry is broken by a tunneling field in a first order 
phase transition, and slow-roll inflation inside the nucleated bubbles is driven by the pseudo-Goldstone boson. In 
section III, we study the quantum fluctuations of the pseudo-Goldstone in the bubble background. Section IV is 
the core of the paper, where we argue that after tunneling, we do not obtain an inflnitc open universe, but an 
infinite ensemble of quasi-open universes inside a single bubble. Section V is devoted to more general models, like the 
"coupled" and "uncoupled" two-field models. In Section VI we briefiy describe the observational implications of our 
results and in Section VII wc summarize our conclusions. 



'^In the case V" <^ the phase transition can proceed via the Hawking-Moss instanton However, this channel represents 
tunneling to the top of the barrier of a region of size , and does not lead to an open universe. 
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FIG. 2. The inflaton potential for the supernatural model. An approximate U(l) symmetry is broken through bubble 
nucleation. We call a the tunneling direction and (j> the direction orthogonal to it. The instanton fo(r) interpolates between 
false vacuum at r ^ oo, and (Tb at r = 0. 



An attractive scenario for open inflation is the model of a complex scalar field with a slightly tilted mexican hat 
potential, see Fig. ^, where the radial component of the field does the tunneling and the pseudo-Goldstone does the 
slow- roll. This model was called "supernatural" inflation in Ref. because the hierarchy between tunneling and 
slow-roll mass scales is protected by the approximate global U(\) symmetry. 

The action is given by 



where Vq is C/(l) invariant and Vi is a small perturbation that breaks this invariance. It is asumed that V has a local 
minimum at $ = which makes the symmetric phase metastable. We shall consider a 'tilt' in the potential of the 
form Vi — A^{(t)G{(I)) where A is a slowly varying function of a which vanishes at it = 0. For definiteness we can take 
G= (1- cos </)/«). 

The idea is that a tunnels from the symmetric phase cr = to the broken phase, landing at a certain value of (f> 
away from the minimum of the tilted bottom. Once in the broken phase, the potential Vi cannot be neglected, and 
the field (f> slowly rolls down to its minimum, driving a second period of inflation inside the bubble. Another attractive 
feature of this model is that depending on the value of (j) on which we land after tunneling, the number of e-foldings 
of inflation will be different. Hence it appears that in principle we can get a different value of the density parameter 
in each nucleated bubble. As we shall see, however, this picture is somewhat oversimplifled. 

We should point out that the supernatural model is not free from certain restrictions. Indeed, in order for the 
pseudo-Goldstone to realize inflation as in the simple free field "chaotic" scenario, we would need v > Mp, where Mp 
is the Planck mass. On the other hand, if Vq is a typical quartic potential, the bubble walls would undergo topological 
infiation for v > Mp, and this would spoil the open scenario. Topological infiation occurs when the thickness of 
the walls is larger than the Hubble rate at the top of the potential barrier separating two local minima (degenerate 
or not). This is the same condition under which the Coleman-de Luccia instantons cease to exist, and we have a 
Hawking-Moss transition instead. Hence the condition v > Mp also represents the regime where the transition is 
not of the Coleman-de Luccia type, as it would be necessary for a successful open universe. As emphasized in Ref. ||^ 
in the case of single-field open inflation, a transition of the Hawking-Moss type would leave unacceptable anisotropics 
in the CMB. These constraints can be made less severe by choosing a suitable form for Vq, with higher curvature at 
the top of the potential between the two minima, or perhaps a special form for G. In any case, we shall take v ^ Mp 
in what follows. 



II. SUPERNATURAL INFLATION 




(1) 



where we use the metric signature (— , +, +). Expanding the field in the form $ 
the expectation value of a in the broken phase, we consider a potential of the form 



{(T/^/2)exp{i(j)/v), where v is 



V^Vo{a)+Vi{a, 0) 
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III. QUANTUM FLUCTUATIONS 



In this section, we compute the amphtude of quantum fluctuations of the pseudo-Goldstone field in the supernatural 
inflation model. For this, we need to review the formalism for quantizing fields in the background of a bubble. 
The field equation for <I> is 

□ <I>-^ = 0, (2) 
where V = Vq-\-V\. In terms of the modulus a and phase , we have 



dV(a, (t>) V dVia, < 

+ i TT— 



(3) 



da a 

It should be noted that classical solutions with (j) — (j) — const, exist only if is an extremum of G, so that 

=0 

This is true also for the Euclidean solutions (instantons) describing the tunneling, so strictly speaking these instantons 
can only take us from the false vacuum to the extrema of Vi in the broken phase. Even if we relax the condition 
that (j) should be constant, there may not be any instantons which can take us to non extremum values. Of course, 
this nonexistence of the corresponding instanton does not mean that the field cannot tunnel to a non-extremum value 
of ^, it simply means that tunneling away from the extremum will be somewhat suppressed. This question will be 
addressed in Section IV. 

Since the effect we are studying is not due to gravity, we shall start with the case of a bubble in flat spacetime. 
Including gravity is quite straightforward and will be done below. According to the theory of vacuum decay [Q, 
the tunneling rate is dominated by the 0(4) symmetric solution of the Euclideanized equations of motion (^) with 
appropriate boundary conditions, which is called the instanton or bounce, and which we shall write as 

(T = cro(T), (p = (p — const, (4) 

Here, we have introduced the Euclidean radial coordinate r = (X^ + T^)^/^, where (Te, X) are cartesian coordinates 
in Euclidean space. As we move from spatial infinity to the origin, the bounce interpolates between false vacuum 
cto(t — !■ oo) = and a certain value of the field in the basin of the true vacuum, see Fig. |^, cro(O) = <Tb. In addition, 
the bounce has to satisfy the boundary condition i7o(0) = 0. The solution describing the bubble after nuclcation is 
given by the analytic continuation of the instanton to Minkowski time T through the substitution Te = —iT. Then, 
the bubble solution depends only on the Lorentz invariant 'distance' to the origin (X^ — T^)^/^, where (T, X) are the 
usual Minkowski coordinates. 

It is useful to change to the new coordinates 

r = (x2-T2)i/2, p^tanh-i(r/|X|), (5) 
in terms of which the line element reads 

ds^ = dr^ + a%{T)dnds^] (6) 

Here dfl^s^ = —dp^ + cosh^ p{d9^ + sin^Od^p^) is the line element of a 2+1 dimensional de Sitter space of unit 'radius', 
and in flat space a^ir) — r. Including gravity, ue has to satisfy the (Euclideanized) Friedmann equation, as described 
below. This 2+1 dimensional de Sitter space can be thought as the hyperboloid swept by the bubble wall during its 
time evolution (suitably rescaled). In spite of its name, the coordinate p is timelike, whereas r is a 'radial' spacelike 
coordinate. 

The above coordinates cover only the exterior of the light-cone from the origin. In order to cover the interior, which 
is where the open universe sits, we use the coordinates 

t = (T2-X2)i/2, r = tanh-^(|X|/r). (7) 

In terms of these the metric reads 

ds^ = -di^ + a^(t)dfijj3, (8) 



4 



where dfi/fs = dr^ + sinh^ r{d9'^ + sin^Odip^) is the metric on the unit 3-diniensional hyperboloid and (in flat space) 
the scale factor is given by a(t) — t. 

Notice that (^) is the metric of an open FRW. When gravity is included, the scale factor a{t) will no longer be 
proportional to the cosmological time t but it will be given by the solution of the Friedmann equation 

In the general case, the metrics and (||) are related by the analytic continuation of coordinates and scale factor in 
the following way 

t — —ir, r — p + i — , a{t) — —iaE{it). (9) 

These relations can be used to analytically continue solutions from the outside to the inside of the light-cone from 
the origin. 

In what follows, we shall assume that tunneling occurs along the real direction for $, i.e. = in (^), and we shall 
consider perturbations around the classical solution of the form 

V2$ = cro(T) +(^1 +i(p2. (10) 

Note that 

<j^^v^. (11) 
Substituting into the action (|l|) we obtain the second order action for linearized perturbations 

= So[ao\ + j d^x^pao ~ V^{(To)]<^i + Si[vi] + 52^2], (12) 

where 

Si^~\ I d^x^g {d^^id^vi + V^'iaoM) , (13) 



S2 



TO^(a-o) 



Here 



m2(ao) = ^A4(ao)G"(0)= 



A^^o) 



(14) 



(15) 



is a small r-dependent 'squared mass' due to the potential V\. In the last equality we have used G of the form 
G — {\ — cos<j)/v). When the field a is in the broken phase, then m is the mass of the pseudo-Goldstone. Of course, 
using the unperturbed equations of motion the linear term in ([12| ) drops out. 

The action ( p^ ) has been studied in some detail in the past [£l[,Q because it is the same as the one for a one-field 
model. In particular, it describes the fluctuations of the bubble wall itself. Here, we shall concentrate on 5*2 [(^2], 
which describes fluctuations in the direction transverse to tunneling. 

In order to study quantum fluctuations, the field ip2 is expanded as a sum over modes times the corresponding 
creation and anihilation operators 



V2 



(16) 



The equation of motion satisfied by the modes is 



□ 



K(^o) 



+ m^((To) 



fplTU — 0. 



Following |2llJl3], we take the ansatz 



^pim = a£;^(T)Fp(r)3^p/,„(a;*), 



(17) 
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FIG. 3. The effective potential in Eq. (^l|). In flat space (dashed line) the potential grows without bound at large rj, whereas 
including gravity it tends to zero. At ry ^ — oo, which corresponds to the center of the bubble, Vcs tends to zero. 



where — {p, 6, ip) are coordinates on the 2+1 de Sitter space spanned by the motion of the bubble waU. Introducing 
the conformal coordinate 77 defined through the relation aE{T)drj — dr, the equation of motion separates into a 
Schrodinger equation 



dr] 



ni ((To) 



n 

6" 



(18) 



where the separation constant plays the role of an energy eigenvalue, and a Klein-Gordon equation for the modes 
of a scalar field of mass + 1 living in the 2+1 dimensional de Sitter space 



(3) 



(19) 



Here '^'^^D is the covariant d'Alembertian in this lower dimensional space, and TZ is the four-dimensional Ricci scalar, 
which from the unperturbed Einstein's equations can be written as 7?. = 87rG'[4T^(ao) + (CTg/as)^]. 

The modes ippim should be Klein-Gordon normalized on a Cauchy surface such as p = 0. This amounts to Klein- 
Gordon normalizing the lower dimensional modes 3^p/m in the 2+1 dimensional sense, and then normalizing Fp as in 
the Schrodinger problem pl|]. 



FpFp'dr] = 5, 



pp' ' 



(20) 



where the delta function will be discrete or continuous depending on whether we are considering discrete or continuous 
eigenstates of the Schrodinger equation. In flat space, ue — Ro&^ , where i?o is an arbitrary constant which can be 
conveniently taken to be of order the radius of the bubble at the time of nucleation. Therefore, the effective potential 
in the Schrodinger equation 



hm (cro) - — 

do D 



(21) 



tends to zero aX rj ^ —00 (center of the bubble) and to infinity at 77 — > +00 (false vacuum), see Fig. ^. In curved space, 
it can be shown that a£;(?7 — > ±00) 0, so VcS vanishes at both ends. Therefore, in both cases, the spectrum will 
be continuous for p^ > 0, and there may be a discrete spectrum for p^ < 0. 



A. Quantum state of a nucleating bubble 



In a time dependent background, the choice of a vacuum state is always somewhat ambiguous 22 1. Here, this 
ambiguity corresponds to the freedom of choosing the 'positive frequency' modes ypim on the hyperboloid. In principle, 
the ambiguity can be resolved dynamically if the initial quantum state before the bubble nucleates is given. 

The quantum state of a nucleating bubble has been extensively studied both in flat and in curved space [ p3[ . In 
our model, (^92 is treated as a free field which couples to the bubble via a a-dependent mass term. For this type of 
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models, it has been shown that if the initial quantum state is de Sitter invariant before the bubble nucleates, then 
right after bubble nucleation the field Lpi will be in an 0(3, 1) symmetric state. This is perhaps not too surprising: 
the appearance of the bubble breaks the 0(4, 1) de Sitter symmetry by selecting a "nucleation point" in spacetime, 
but otherwise the bubble solution respects an 0(3, 1) subgroup of isometrics. 

What this means is that the positive frequency modes 3^p/m must be taken as the Bunch-Davies modes, which 
guarantee the desired symmetry. These are given by 



r(z + 1 - iv)V{i + 1 + ip) 



1/2 p 



[i smh. p) 



\J i cosh p 



Yin 



(22) 



where P/,' are the Legendre functions and are the usual spherical harmonics. When analytically continued to the 
inside of the light-cone, through the relations (0), they become 



V sinh 1 



-Y, 



Ira • 



(23) 



These are proportional to the often used harmonics Ypim which are normalized on the hyperboloid H 



_r(zp) r([+i_^ 

^/2 Vr(^ + 1 + *P) 

These analytically continued modes are normalizable on only for > 0. Since the mode with p^ = has wave- 
length comparable to the curvature scale, the non-normalizable modes withp^ < have been dubbed "supercurvature" 
modes |9|. Writing p = —iA, the supercurvature modes are given by 



r(A + / + i)r(-A + z + i) 



1/2 p 



-1-1/2 



A-1/2 



(cosh r) 



V sinh? 



(24) 



(We have added a comma after the subindex A to indicate that it is the value of ip rather than p). We repeat, 
however that if the corresponding bound state exists in the Schrodinger equation (|l^), then these modes are perfectly 
normalizable on the Cauchy surface, and hence they must be included in the expansion of the field operator (16). 



B. Degenerate case 



To begin with, let us neglect the mass term ([T^) which comes from the tilt in the potential, Vi 0. Then, using 
the equation of motion satisfied by ctq , it is straightforward to show that 

Fi = Afasao (25) 

is a solution of dlq) with eigenvalue = — 1 , or A = 1 . Moreover, this solution is normalizable and it belongs to the 



discrete spectrum. The normalization constant M is found from ( ^ ) 



aE(rl{T)dT 



-1/2 



(26) 



Here, we have changed back to the physical coordinate r, which measures the physical distance to the center of the 
bubble. In flat space Tmax is actually infinite, but the integral is finite because (Tq vanishes exponentially fast outside 
the bubble. Including gravity Tmax becomes a finite value, so the integral is also finite. The explicit value of M can 
be calculated numerically for any given model. If i?o is the size of the bubble and (Jb is the value of the field at the 
center of the bubble, see Fig. |l|, at the time of nucleation, then we can estimate 



V2 
Roab 



Tunneling rates in the case when there is an exact internal symmetry have been recently investigated in |24| 
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The above estimate is for bubbles such that Rq is small compared with the Hubble radius, so that qe ~ t. Also, 
we have substituted (Tq by Ub in the integrand, and we have integrated r from zero up to Rq. Because (Tq < o";, our 
estimate is actually a lower bound on Af, and in some thick wall models the effect can be somewhat higher. 

If we asume a large energy difference between false and true vacua, then we are in the thick wall regime and Rq is 
of the order of the thickness of the bubble wall 

Rq - M~\ 



Here M is the mass of the field in the false vacuum. If, on the contrary, the vacua are sufficiently degenerate, then 
we are in the thin wall regime, and the expression for Rq can be found e.g. in Appendix C and Ref. [Eq]. 
The normalized mode has an amplitude 



'Pl,lm J^l.iji 



.{Xi). 



(27) 

, but a distribution 
the X* dependence 

in ( |27|) shows that, after nucleation, different points on the bubble have different values of the angle (j). 

When analytically continued to the interior of the light-cone from the origin, through the relations t is replaced 
with the cosmological time and <7o(t) quickly follows its evolution towards its expectation value ctq ^ v > ai,. Hence, 
inside the light-cone, the normalized supercurvature mode will take the form 



Physically, what happens is that the field does not simply tunnel to a sharply defined value of q 
of values. Taking into account that the phase of our complex scalar field $ is given by Eq. (nw 



Hq at 



(28) 



which is time independent. In order for the bubble solution to exist, the mass of the field in the false vacuum should 
be M ^ Hp, where Hp is the Hubble rate in the false vacuum. In our model. Hp is considered to be much greater 
than the Hubble rate in the true vacuum ^^t-]^ Therefore, the amplitude of supercurvature perturbations is of order 
(v/<Th)M . This exceeds by far the amplitude of the usual "subcurvature" fluctuations, which is only of order Ht- The 
corresponding effect in the CMB places severe constraints on the model We shall come back to this question in 
Section VI. 

Note that for A = 1 the amplitude of the homogeneous mode I — diverges because one of the gamma functions 
in (|^ has vanishing argument. This is related to the fact that, strictly speaking, the I = mode should have been 
quantized as a collective coordinate instead of as a harmonic oscillator ||26| . The divergence simply means that all 
values of cj) are equally probable after nucleation. When we include the effect of Vi, the degeneracy will be broken 
and the I = mode will have a finite amplitude. 



C. Non-degenerate case 

When the tilt Vi in the potential is included, = — 1 is no longer an eigenvalue of the Schrodinger operator (p^). 
However, it is clear that for small there will still be a discrete eigenmode whose eigenvalue we can calculate in 
perturbation theory. Denoting by | — 1) the unperturbed bound state (Pq), the perturbation to the eigenvalue will be 
given by 

7 = + 1 = (-l|a|m2| - 1) = y a%alm'^{ao)dr]. (29) 

Again, 7 can be computed numerically for any particular model, but we can estimate it as being of order 

/a|a^m^(go)rfT 1 ^2 /on^, 
j aEcr^[T)dT 2 



^ We are considering here strongly non-degenerate minima. However, in Ref. g they also consider various depths of the 
central minimum, depending on radiative corrections, and in some cases (5* = 327r^A) the two minima become degenerate, 
Ht ~ Hp- This weakens the constraints and makes the model viable in certain range of parameters. 
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where = A^(crb)/cr^, see Eq. (|l5|). The above estimate uses the same approximations as the estimate for M after 
Eq. (H). 

The normaUzed I — mode will now be given by 

7 ^ Ar ^1 Mt) 1 sinh[(l - 7)V^r] 

tt/cq iJb 7 ' smhr 

where, from now on, wc shall use the notation ip'^ to denote the modes with = — 1 + 7. In this expression, and 
the ones that will follow, we give only the unperturbed time dependence, without including the first order correction 
in 7. The field ao(t) quickly rolls down from ab and settles down to its minimum at cto = v. It is understood that, 
after that, the time derivative of the mode will be dominated by the corrections of order 7, which slowly drive the 
Goldstone modes to the minimum of the tilted potential. The amplitude of the higher-Z excitations will also suffer 
corrections of order 7, but since these amplitudes were already finite, the effect on those is not so dramatic. To leading 
order, the amplitudes are the ones calculated in the previous subsection. 

It is interesting to look at the distribution of the field on a hyperboloid t = const. Note that the amplitude of the 
I = mode near the origin r = is of order 

which is a factor 7"^/^ larger than the amplitude of the individual I > 1 modes found in Eq. (p^). But the amplitude 
of the I = mode decays exponentially for r 3> 7~^, which means that at large distances it will become negligible. 
However, the quantum state that we have chosen is 0(3, 1) symmetric, which means that the r.m.s. fluctuation of 
the field cannot depend on r. Therefore, the loss in amplitude of the I = mode as we move away from the origin 
has to be made up for by the joint contribution of the I > modes, smeared over a suitable length scale. This is 
analogous to what happens for a massive field in de Sitter space, except that here we are considering a spacelike 
manifold t = const, rather than a spacetime. In the Appendix B we show that for r ^ we have 



1=1 ^ ^ 



where Z, 3> 1 is a certain cut-off. If we smear over a fixed comoving length ^, as we move away from the origin we 
have to include more and more modes in the sum. Since the l-th multipole has wavelength proportional to (sinhr)/Z, 
we take Z, = sinhr/^. With this, we find {Aip'^)'^ « (A/'^CTQ/27r^7)(2^)^^. Notice that the result is rather insensitive 
to the choice of ^. As long as | ln^| ^ 7~^, the added contribution of all relevant modes at large r is the same as 
the contribution of the I = mode near the origin, given by (^) . We also show in Appendix A that the two-point 
correlations on a, t = const surface die off with comoving distance d as e~^'^l'^ . 

Hence, around the time when the field ctq settles down to its minimum ^ v the supercurvature fluctuations of 
the field are of order 

[^V )~ ^ ^^J ^^^2 ^ A2(afc)mi?2 ~ rn ' ^ ^ 

Here, m = K^{v)/v is the mass of the pseudo-Goldstone in the true vacuum. In the last step we have asumed that 
A(cr) is a slowly varying function of cr, and set i?o ^ . Note that the fluctuation in if 2 can easily reach Planckian 
values. It suffices to take M ^ 10^^ GeV and to ~ 10^^ GeV. With these values, the field is displaced enough from 
its minimum that it can drive inflation (recall that we are asuming v ~ Mp). 



IV. QUASI-OPEN UNIVERSE 



Tunneling to a large value of the field is usually understood in an "adiabatic" sense. The idea is that since the 
motion of the phase is dictaded by the explicit symmetry breaking potential Vi, it will be much slower than the 
motion of the radial component of the field dictated by the large U(l) symmetric part of the potential. Hence, one 
can estimate the rate for tunneling at any (f) by solving the Euclidean equations of motion for a while </> is kept 
as a frozen parameter. This frozen parameter is then used as the initial value of the slow-roll field on the t = t^, 
hypersurface. Here, as above, t^, is the time at which ag reaches its expectation value v. 

However, it is clear that when the spread of the pseudo-Goldstone on the t = t^, hyperboloid is comparable to its 
range, {Af>'') > nv, the picture that each bubble nucleates with a different value of (p2 is not adequate. Instead, all 
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FIG. 4. On at — const, hypersurface inside the bubble, the co-moving coherence length of the slow-roll field is r --^ 7~^- If 
the r.m.s. fluctuations of the field are large, then regions where the field is large and positive will be next to regions where the 
field is large and negative. These will be separated by regions where the field is small and there is no second stage of inflation. 



of the 'vacuum' manifold is pretty much sampled inside a single bubble. In this case, inflation does not take place 
coherently on the t = const, hyperboloids. As mentioned in the previous Section, the comoving correlation length for 
quantum fluctuations (in units of the curvature scale) is d ~ 7"^. Thus, patches of comoving size d ~ 7"^ where (p2 
is large and positive would be right next to patches where ip2 is large and negative. These patches would be separated 
by regions where ip2 is small and the universe does not inflate. Patches with positive and negative values can also 
be separated by "walls" where (p2 is close to ttv. In the case v « Mp these domain walls of the pseudo-Goldstone 
potential would be topologically inflating |^^. Hence, instead of a smooth inflating hyperboloid, we have a patchy 
mosaic of inflating regions, as depicted in Fig. ^. In principle, each patch can give rise to a successful cosmology, but 
this cosmology will not be an open universe in the traditional sense. At best it will be a 'quasi-open' one, i.e. one 
which locally resembles an open FRW. 

When (A(/3^) <C nv the spread in the distribution of the pseudo-Goldstone in a t = section is small. What this 
means is that if the nucleated bubble is described by the 0(3, 1) symmetric quantum state, then most of the surface 
t = t^ has a non-inflating value of the field. However, in an inflnite hypersurface, there will be a certain density of 
occasional large fluctuations which will lead to inflating islands of comoving size d ^ 7^^. Each one of these islands 
will be a quasi-open universe. 

This is in clear opposition with the conventional adiabatic picture described in the first paragraph. It is interesting 
to pursue the adiabatic picture for a moment, in order to see in which sense it is adequate or not. To keep the 
discussion simple, let us consider tunneling to a range of values of (f) for which the linearized expressions (^2|) are still 
valid, but sufficiently large that it can be distinguished from tunneling to the bottom = 0. This will be the case if 
ip2, see Eq. (pT|), is in the range 



(A(y92) < 1^2 < 'f , 



(34) 



where (Aiy92)^ = {N"^ 0^/211^^) was computed in the previous section. 

Since ipi and ip2 decouple, we may take an approximate Euclidean solution of the form $ = (1 + i(j)/v)ao{T)/\^, 
where is taken as constant in the adiabatic approximation. Substituting this configuration in the Euclideanized 
version of ( p^ ) we find, after straightforward algebra. 



Se[^] 



77r 



If the decay rate is proportional to e "^^ , the relative probability of having a bubble with a certain value of (f2 at 
nucleation will be 



P(V?2) ~ exp 
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2{Aip2y 



(35) 



where we have used ( [ll| ) and, perhaps not too surprisingly, found {Aip2)^ — (A/'^crQ/27r^7), the same expression 
obtained in the previous section from considerations of quantum fluctuations in the 0(3, 1) invariant state. Thus, 
two approaches which in principle are aimed at answering different questions end up giving the same answer. Here, 
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we were asking how likely is it for a bubble to nucleate at a large value of 0, whereas in the previous Section we were 
computing the amplitude of fluctuations inside a given bubble. Undoubtedly, both questions are related, since at least 
locally we cannot distinguish a large value of the field induced by the nucleation of the bubble from a fluctuation of 
the field inside the bubble. 

However, even though the adiabatic approximation may give the right answer for the probability of tunneling to a 
large value of 0, it suggests the wrong picture for bubble nucleation. Since = const is used in the above estimate, 
one might imagine that an inflnite open universe with homogeneous Lp2 in the range (^) can be created. We shall 
argue that the probability for this to happen is actually zero. First of all, the formula P ~ exp[— S's] is only justifled 
when the Euclidean action is evaluated on a solution of the equations of motion. But cj) — const is only a solution 
in the case when Vi is neglected. We can try and correct this configuration so that it will be a solution, while still 
keeping 0(4) symmetry. In the linear regime, what this means is that we want a solution of Eq. ( p^ ) with = — 1, so 
that Eq. (M) is satisfled by ypim{x^) = const, and we have a homogeneous solution inside the bubble. However, for 
> the lowest eigenvalue isp2 = — 1+7>— 1, which means that the solution with = —1 is not normalizable. 
As a result, the corresponding Euclidean action is badly divergent. If the action is regularized with a cut-off and if 
we take the = —1 solution to be well behaved at the center of the bubble t — 0, then it is easy to show that the 
action starts growing exponentially, Se ~ exp[2M(Tc — Rq)], as the cut-off Tc in the radial direction r becomes larger 
than the size of the bubble. Here M is the mass of $ outside the bubble. Hence it is not justified to say that the 
estimate ( |35| ) gives the nucleation rate for a homogeneous bubble. Rather, using a homogeneous solution we would 
get a divergent action and hence a vanishing probablility. 



A. Creation of a quasi-open universe 



To compare, we can now ask what is the amplitude for tunneling from false vacuum to a spherically symmetric but 
inhomogeneous configuration with a large value of ip2 inside the bubble. This is what we call an "inflating island" or 
quasi-open universe. Again, this amplitude will depend on the action of a semiclassical Euclidean trajectory. In order 
to make the metric into one of Euclidean signature, we must consider the analytic continuation of the coordinates 
t = —IT and r = ite- With this we have 

ds^ ^ dr^ +a|(T)drj53, (36) 

where dfi^s = dr^^ + sin^rEidO^ -\- sin^Odfj)^) is the metric on the three-sphere, and the range of te is from to it. 

The semiclassical trajectory we shall consider is simply the analytic continuation of the / = supercurvature mode 
0® 

V2^A f{T)g{rE), (37) 

where /(r) — ao{T) + 0{j) and 

sin[(l — ^y^^rE] 



gi^E) 



smr^. 



Note that this solution is not regular at one of the poles of the three-sphere, rE — tt. This is of course expected, 
since g is an eigenfunction of the Laplacian with eigenvalue 7, which is not in the spectrum. Hence, ( p7| ) is not a 
regular "bounce" around which we should expand in order to find the false vacuum persistence amplitude This 
is not a problem for us here, because we are not estimating the total decay rate of the false vacuum but a particular 
transition amplitude. For this, we only need "half" of the Euclidean solution, which interpolates between false vacuum 
at T — > 00, rE ~ and an inhomogeneous configuration at rE = 7r/2. At rE = 7r/2 this solution is matched to a 
Lorentzian solution at p ^ (the spacelike surface where the bubble nucleation takes place) simply by the analytic 



"^The semiclassical trajectory we consider does not have vanishing "temporal" derivative at the bounce point rs = 7r/2 where 
we match the Euclidean solution to the Lorentzian one. Hence it is a complex trajectory which has a small imaginary part 
(of order 7) in the classically forbidden region. The imaginary part decreases exponentially fast in the Lorentzian section on a 
timescale of order Rq^ 

^This is a blessing, because an inhomogeneous instanton would cause the well known problem that the decay rate should be 
multiplied by the infinite volume of the Lorentz group. 
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continuation t = —ir and p = i{rE — 7r/2) [see Eq. (^]. The solution is then propagated to the interior of the bubble 
(the "open" universe) through Eq. @. 

The transition amplitude is WKB suppressed only in the Euclidean regime; the Lorentzian evolution contributing 
an oscillatory phase. Therefore 



where 



Sk = 27r 



dT a%{T) 



r/2 



dr 



E sm ve 



1. (^\ 

«l \drEj 



CTO 



It is easy to show that since ip2 is a solution of the equations of motion, after integration by parts the only contribution 
to the integral comes from the boundary term atr£; = 7r/2. A straightforward calculation gives 



Se^A^'-^ 



J alaldi^^A^ 



77r 

W2 



where A is de constant introduced in (|37|) . After analytic continuation to the interior of the bubble, we have tf2 {r ^ 



by 



Aao. Hence, the relative probability for nucleating through the inhomogeneous trajectory (pTh is again gi" 



1*1 



exp 



2 

-V2 



(38) 



Although this is in perfect agreement with ( pq ) , it is now clear that it doesn't mean the whole open universe will have 
the value ip2, but only a patch of comoving size ~ will have this value. 

Thus, tunneling to a value of the field which is far from the one indicated by the 0(4) symmetric instanton is 
perfectly possible, with a somewhat suppressed probability. However, the resulting universe is not an infinite open 
universe but just a quasi-open one. 



B. Many universes in one bubble 



The arguments used in the prcvios subsection leading to Eq. ( p8D are somewhat heuristic. In particular, we have 
not attempted to justify why the semiclassical trajectories of the form (^^ should be the only relevant ones. Note, 
however, that the probability distribution (38) for nucleating at a high value of the field ip2 near r = 0, is the same as 
the Gaussian distribution for the amplitude of the ^ = 0, = — 1 + 7 mode in the 0(3, 1) invariant state. In fact, the 
possibility of nucleating at different values of the field is already accounted for by this quantum state and need not 
be considered separately. The analysis of Refs. ||2^], whose result we described in Section HI. A, takes into account all 
paths in t^-field space, not just the semiclassical one used above. That analysis should be regarded as a more rigorous 
derivation of the result ( pSj ) . 

Because of the invariance of the quantum state, it is clear that there is nothing particular about the point r — 0, 
and an inflating 'blob' is equally likely to develop around any point on the t — t^, hyperboloid. Therefore we are 
led to the picture where in each comoving volume of size comparable to the correlation lenght 7~^, the probability 
distribution for ip2 is given by (|38|). Clearly, since the volume of the hyperboloid is infinite, inflating islands with all 
possible values of the field at their center will be realized inside of a single bubble. We may happen to live in one of 
those patches of comoving size d ^ 7""'^, where the universe appears to be open. 

Also, in the supernatural model, it is possible to modify the shape of the potential near the false vacuum so that 
there is a misalignment between the prefered direction for tunneling Q and the direction of the minimum of the 
pseudo-Goldstone potential in the broken phase. In this case, we expect that the t — t^, surfaces inside the bubble 
will have a mean value of = 0c 7^ 0, determined by the most probable escape path (i.e., the instanton, which in 
this case will not land on = 0). This value of 0c will determine the number of e-foldings of inflation and hence 
the mean value of the density parameter fig on the hyperboloid. Let us call this value flc- If the tunneling path is 
not too narrow, there will still be a supercurvature mode which will cause fluctuations in the density parameter on 
co-moving scales of order 7"^, which are of course much larger than the Hubble radius. The picture is then that 
we have an ensemble of large patches with different values of the density parameter. This is an interesting situation 
which deserves further study. However, since this model involves more parameters, for the remainder of this paper 
we shall concentrate on the simplest case discussed above, where the preferred tunneling direction and the minimum 
of the pseudo-Goldstone potential are aligned. 
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V. MORE GENERAL MODELS 



In this section we shall consider the class of two-field models with a potential of the form 

Models of this type were introduced in Ref. Here Vb is a non-degenerate double well potential, with a false vacuum 
at cr = and a true vacuum a,t a — v. When a is in the false vacuum, Vb dominates the energy density and we have 
an initial de Sitter phase with expansion rate given by « (87rG/3)Vo(0). Once a bubble of true vacuum a = v 
forms, the energy density of the slow-roll field cf) may drive a second period of inflation. 

As pointed out in Rcfs. ||^, the simplest two-field model of open inflation, given by ( p9| ) with g — Q and m ^ 0, is 
actually a quasi-open one. Since there is no coupling between the two fields except for the gravitational one, we shall 
call this the "decoupled" two-field model. In this model, infiation starts chaotically at large values of V{a, 0) ^ M^. 
In some regions of the universe, the field a will be trapped in the false vacuum, while </> rolls down from large values. 
If a bubble nucleates at a point where (/) ~ Mp, the value of the slow-roll field will be large enough to drive a short 
second period of infiation inside the bubble. One problem with this model, is that the slow-roll field moves also 
outside the bubble, so the synchronization of the = const, and a = const, surfaces inside the bubble is not perfect, 
as pointed out in Ref. 

In Ref. H, the classical evolution of the slow-roll field from the outside to the inside of the bubble was studied, 
and it was found that on the hypersurface t ~ Hp^ the high value (j) ~ Mp decays exponentially with the distance to 
the origin as 

4> cx exp[-7cr/2] (40) 

where 7c « 2mjp/3Hp. Hence, the size of the infiating region in this model is finite. We use the subindex c to stress 
that this result follows from purely classical evolution. The larger Hp, the larger will be the infiating region. The 
reason is that the cosmological friction term in the equation of motion for cf) is proportional to Hp, so the larger is Hp, 
the slower the field (j) will roll down the potential outside the bubble, and the better is the synchronization between 
(j) = const, and cr = const, surfaces. However, Hp cannot be taken to be too large because otherwise the quantum 
fluctuations generated outside the bubble produce too large an amplitude for the supercurvature mode inside the 
bubble H] . The combination of these two effects severely constrains this model ||l^ . 

In order to construct a truly open model, Linde and Mezhlumian suggested taking m — and g 0. We shall 
call this the "coupled" two-field model. In this way, the mass of the slow-roll field vanishes in the false vacuum, 
and it would appear that the problem of classical evolution outside the bubble is circumvented. However, this is not 
exactly so, and the whole class of models ( |3^ leads to quasi-open universes. The basic reason is that, as we shall see 
below, the (linear) equation of motion for (j) in the presence of the bubble, does not admit 0(3, 1) invariant solutions 
which are regular at the origin, except for the trivial one, = 0. Thus, we are back to a situation analogous to the 
supernatural infiation model. 

Even if the mass of the field in the false vacuum vanishes, one must not expect that will not evolve at all outside 
the bubble. Fig. || shows the result of a numerical evolution of the field (j) in the coupled model. The figure represents 
a conformal diagram of a bubble expanding in de Sitter space (for simplicity, the gravitational field of the bubble 
has been neglected). The bubble wall is indicated by the thick timelikc hyperbola. As initial conditions, we have 
taken (j) = const, and (f> = 0. Surfaces of constant are indicated by different shadings. Even though the field is 
massless outside the bubble, we find that it does not stay exactly constant there. Due to the finite size of the bubble 
at the time of nucleation, the field cj) feels the presence of the bubble everywhere inside the light-cone from the "point" 
P. As a result, inside the bubble, the hypersurfaces (f) — const, are not perfectly synchronized with the a = const. 
hypersurfaces. Thus, we are back to a situation where the infiating region inside the bubble has a finite size, as in 
the decoupled model (note that the = const, lines cross the bubble wall trajectory). 

Note that this effect is due to the finite size of the bubble. We shall see below that the effect is of order [HpR^Y. 
In Fig. ^, the parameters have been chosen so that the effect is very dramatic and the size of the inflating islands is 
comparable to the curvature scale, but one can choose parameters so that the inflating islands are as large as desired. 
However, except in the case where gravity is neglected so that Hp = 0, their size is always finite and the large value 
of the field at the time of nucleation ends up decaying at large distances r from the origin. 

Hence, just as in the case of the supernatural model, the infinite t = const, surfaces would be almost empty at large 
distances, if it wasn't for the occasional quantum fiuctuations which may ignite inflating islands here and there. 
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FIG. 5. Conformal diagram of de Sitter space with a bubble expanding in it. Here rj and p are the usual conformal 
coordinates in the closed chart. The bubble wall is represented by the curved thick line starting at the point P. The figure 
shows the result of a numerical evolution of the slow- roll field in the coupled model ( |39[ ) with mp = and mr ~ 1.5Ht- The 
self-gravity of the bubble has been ignored and we have taken Ht ~ Hp. The initial conditions for at 77 = 7r/2 are (j) = const 
and = 0. In spite of the fact that the field is massless outside the bubble, it starts evolving everywhere inside the light-cone 
from the point P. As a result, the surfaces of constant 6, separating regions with different shadings, are not well synchronized 
with the t = const, surfaces inside the bubble (see Fig. hi). In the plot, the field decays by one tenth of its initial value between 
consecutive (f) = const, lines. 



A. Quantum fluctuations 



As in the case of the supernatural model, we expand the field operator 
operators, 

= ^ fphnaplm + h.C. 

In the present case, the equation of motion for the modes is given by 

OiPpirn - [m^ + gcrQ]ippi„i = 0. 
Here ctq is the "background" bubble solution. Using the ansatz (|l7|), 

we have the following Schrodinger equation for Fp, 



in terms of creation and anihilation 



(41) 



9^1 



n 

6" 



Fp = p^Fp. 



(42) 



This equation determines the spectrum of allowed eigenvalues p^, which correspond to normalizable eigenfunctions 
Fp. All of these eigenvalues have to be included in the expansion of the field operator (41). As before, the harmonics 
ypim must satisfy equation (p^. 

In the case of supernatural inflation, we saw that there was a discrete eigenstate with p^ < 0, which actually 
dominated the r.m.s. fluctuations of the field on a. t — const, hypersurface. We shall see that a similar situation 
happens in this case. 

Since the Hubble rate inside the bubble Ht is smaller than the Hubble rate outside. Hp, we need 



(43) 



in order to have slow-roll inflation inside the bubble. This condition suggests taking a perturbative approach. To 
lowest order, we can neglect the mass term for and the gravitational backreaction of the bubble, so that (p is just a 
massless field in de Sitter space. In this case, Eq. (^) has the well known supercurvature mode with p^ = —1, which 
corresponds to [pi 
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F_i = ^MV) (44) 

The field configurations corresponding to this p are (j) = const, x 3^i,;m- As mentioned before, for I — the harmonic 
3^1,00 is not Klein-Gordon normalizable in the 2+1 dimensional sense. This is because, for a massless field in de Sitter 
space, the zero mode corresponding to the translations (f) (f) + const has to be treated as a collective coordinate and 
not as an oscillator. When the mass of (f> is included, the mode becomes normalizable. Indeed, the bound state will 
shift to a perturbed eigenvalue = — 1 + 7 which, as before, can be calculated perturbatively: 



_ 2 I 1 I 4 



drj. (45) 



The normalized modes corresponding to the discrete eigenvalue take the form LpJ^_^ w {H f / ^/^)ypim{x^) , and in 
particular, for the I — mode. 

The uncertainty of order 7 comes from the fact that we have not evaluated the correction to the "wave function" F_i 
which gives the temporal dependence of the field. This can be done in principle, but it is not really necessary for 
our purposes. It is clear that the mass term will cause the field to have the temporal dependence corresponding to 
slow-roll inside the bubble. 

Near the origin r = the r.m.s. fluctuation of the field A(/) will be dominated by the mode (|46|), and for t — t^ ^ Hp^ 
it will be given by 

A^-f^^^- (47) 

As mentioned in Section III, because of the 0(3, 1) invariance of the quantum state, this will also be the r.m.s. 
fluctuation of the field at any point on the t — hypersurface (see Appendix B). 

In the case of thin walls, the value of 7 can be calculated explicitly (see Appendix C), 

l = \^^+\HlRt{ml~my). (48) 

Here nri^ — w? + gv'^ is the mass of the slow-roll field in the true vacuum, vn?p — nn? is the mass in the false vacuum. 
Hp is the bubble rate in the false vacuum, and i?o is the intrinsic radius of the bubble at the time of nucleation. 

The origin of the different terms in ( ^8| ) is easy to understand. The first one is independent of the existence 
of the bubble, and comes from the fact that the slow-roll field has a mass m?p in the false vacuum (the in the 
denominator can be understood from simple dimensional considerations). The second term is due to the perturbations 
of the effective potential in the Schrodinger equation (|4|) caused by the bubble solution. In the bubble, the scale 
factor aE is of order _Roi so the factor in the integrand of ( |45| ) will yield the factor Rq in front of the second term 

of il). 



B. Inflating islands 

Even though the decoupled model {g = 0) is not a very good candidate to an open cosmological model it is 
instructive to consider it as a first step. In this model, inflation inside the bubble can be initiated because of the large 
"classical" value of </> at the time of nucleation, but at very large distances from the origin r ^ 7^^, the classical field 
dies off, and only the quantum fluctuations remain. 

Let us consider the amplitude of quantum fluctuations. From ( ^ ) we obtain {inp = nix) 

_ 2m| 

Accordingly, from (^7|), the r.m.s. quantum fluctuations of the held on the t = t<t surfaces will be of order 

(A^)«yif— ' (49) 

V 2 27r rnp 
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and the comoving correlation length will be of order 7"^ = the same as we found in the classical case ( ^0| ) from 
a completely different approach. 

In the decoupled model is the same as uit, and Hp is not too much larger than Ht |^,||,|9|. S ince H'p ^ ttij^ , 
we have (Ai/)) <^ Mp and quantum fluctuations on that surface would typically not reach "inflating" values of order 
Mp. Still, as in the supernatural model, an occasional large quantum fluctuation can initiate inflation on a patch of 
size 7^"'". 

Let us reflect upon the meaning of Eq. (^). This r.m.s. amplitude is actually the same as the Bunch-Davies 
one for fluctuations of the slow-roll field outside the bubble ||2^ . Hence, the correct interpretation of the result ( |49| ) 
seems to be the following. Inflation inside and near the bubble wall may start because the field is large at the point 
where the bubble nucleates. However, even after bubble nucleation, the field will continue its random walk outside 
the bubble, and it may occasionally become large. If the bubble wall hits a patch where the field is large, then this 
will generate a local inflating patch inside the bubble, and we might inhabit one of those inflating patches. However, 
this model is not a very good open model, and it would only agree with observations if CIq were very close to one. 

Let us now consider the "coupled" model, where the fields are coupled but mp = ni = 0. In this case 

7 « ^rn'rHlK 

Clearly, by choosing parameters such that the size of the bubble is much smaller than the Hubble rate outside the 
bubble, or such that the mass of the field in the true vacuum is sufficiently small, 7 can be made as small as desired. 
Hence the size of inflating regions can be made as large as desired. In this case, the field is massless outside the bubble 
and quantum fluctuations of pile up to arbitrarily large values far from the bubble. However, from Eq. (|47|), we 
find a finite answer for the fluctuations inside the bubble 

(A</.)«^^-^. (50) 

TT IJItRq 

The first interesting thing to note about this result is that it does not depend explicitly on the Hubble rate inside 
or outside the bubble (the only dependence is through Rq). The second observation is that it is very similar to 
the expression ( p3| ) we had for the supernatural case, so a connection between the physics of both models can be 
anticipated. The finiteness of (^0|) is not surprising, since the slow-roll field is coupled to the bubble, and piling of 
modes in the vicinity of the bubble is suppressed by the mass term. Also, nucleation of bubbles at high values of 
(j) is suppressed because the degeneracy between true and false vacuum is lower. As we discussed in Section IV, the 
quantum state already encodes the information that tunneling to a large value of the field is suppressed. 

A difference with the supernatural inflation case is that now the amplitude of the supercurvature modes with I 7^ 
is of order Hp rather than 1/Ro, hence the constraints on this model from microwave background anisotropics will 
be easier to accomodate. 



VI. OBSERVATIONAL CONSEQUENCES 



Our results from the previos sections have important consequences for two-field models of open inflation. First of 
all, our models are quasi-open, rather than open, which leads to classical anisotropics [|l9|. Second, we saw that in 
the supernatural model, the amplitude of supercurvature excitations is quite large. In this section, we shall give order 
of magnitude estimates for the expected CMB anisotropics from these effects. A detailed investigation of the power 
spectrum will be presented elsewhere [ p8[ . 



A. Classical anisotropies 

Quasi-open universes are finite, and hence they look anisotropic to a typical observer. This effect was studied in 
Ref. for the uncoupled model, and was called a "classical anisotropy" . The name was given because the finiteness 
of inflating islands was due to the classical motion of the slow-roll field outside of the bubble. Clearly, the same effect 
arises in all quasi-open universes we have considered. In some cases the appearance of the inflating island is better 
described as a semiclassical effect, but the resulting inflating islands are just as classical here as they were in Ref. [ p^ . 
Hence, we shall use the same name for this type of anisotropies. 

To proceed, it will be important to distinguish between two different cases. The first case arises when the r.m.s. 
fluctuation of the slow- roll field cj) on the spacelike surfaces t ^ t^, ^ Hp^ is small compared with Mp. In this 
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case, "high peaks" where the field is comparable to Mp and which will lead to inflating regions of co-moving size 
r ~ ^ 1 will be very "rare" on that hypersurface. Here, 7 is the correction to the supercurvature eigenvalue, 
calculated in Section lll.C for the supernatural model and in Section V.A for the more general models. High peaks 
of a homogeneous Gaussian random field tend to be spherical, and so our inflating islands will have approximate 
spherical symmetry. In the opposite case, when the r.m.s. fluctuation of the slow-roll field is comparable to Afp, we 
will have a patchy mosaic of "overlapping" inflating islands, as described in the second paragraph of Section IV. It is 
easy to check that in the second case the "classical" effect is small compared with the effect of quantum fluctuations 
which we shall consider in the next subsection, so here we shall only consider the first case. 

Let us begin with the supernatural model. The quantum state we are considering leads to a Gaussian distribution 
for the random field (j) that is 0(3, 1) symmetric. Hence, to compute the probabilities for the field distribution around 
any point it suffices to study them around the origin r — Q. Here we shall only be concerned with fluctuations due 
to the supercurvature modes, which have a long range. The effect of subcurvature modes can be incorporated in the 
usual way ^. Note that the r.m.s. amplitude for the I — supercurvature mode is a factor larger than the 

amplitude for / > modes (recall that 7 <C 1). Hence, even if the r.m.s. of is far below Afp, there is a certain 
probability for (j) to reach Mp in a certain region near the origin. The spherically symmetric mode is the one that is 
most likely to contribute to this possibility. Even though there is a small probability for this to happen, it is clear 
that only those rare regions with ^ Mp will undergo a second stage of inflation; so they will be the only ones that 
matter. The value of the field on those inflating islands will have the radial dependence of the I — mode, which 
decays as exp(— 7r/2) at large distances, r ^ 1. 

Let us now discuss the more general models where the slow-roll field has a small mass or it is massless outside the 
bubble. In this case, one may ask what happens when a bubble nucleates in a place where the slow-roll field already 
had a large (classical) value. This may occur, for instance, if the whole universe was created at a large value of 4>i and 
at the time when the bubble nucleates (j) is still rolling down from large values. This possibility would in principle be 
relevant for bubbles nucleated at early times, and is the one considered in p9[ |. However, as time goes by, the initially 
large classical value of the slow-roll field in the false vacuum will decrease, and all that will remain are the quantum 
fiuctuations which should be well described by the 0(4, 1) or de Sitter-invariant quantum state. 

Occasionally, fluctuations of the slow-roll field in the false vacuum may create a localized region with a higher 
value of the field. The nuclcation of a bubble on top of one of these regions will not be very different from the case 
discussed in the previous paragraph. Whether the bubble nucleates on one of these high peaks or not, the field outside 
the bubble will continue to fluctuate, and the bubble walls will from time to time bump into regions with a higher 
value of the field, as discussed in Section V.B. Hence, also in this case, there will be an ensemble of inflating regions 
with some distribution inside the bubble. From a formal point of view, notice that the appearance of the bubble has 
selected a point in spacetime, thus breaking the 0(4, 1) invariance, but otherwise respects a residual 0(3, 1) symmetry. 
Therefore it seems reasonable to expect that, at least in a statistical sense, the field inside the bubble will be well 
described by the 0(3, 1) invariant Gaussian distribution, corresponding to the quantum state we have studied. Just 
as in the case of the supernatural model, here we also expect that the high peaks which lead to inflating islands will 
have spherical symmetry, and the value of the fleld on those islands will have the radial dependence of the / = mode, 
which decays as exp(— 7r/2) at large distances from the center of the island. 

The co-moving size of the inflating islands is r ~ 7^^ 3> 1. Since the volume on the hyperboloid grows exponentially 
with the distance to the origin, as sinh^ rdr, a typical observer in a quasi-open universe is most likely to be at r 3> 1, 
where the scalar field behaves as </> oc exp(— 7r'/2). Up to exponentially small corrections, this is the same radial 
dependence that was considered in ||l^. In that case, the fields were uncoupled and 7 — 2mp/SHp. The arguments 
used in to estimate the temperature anisotropics measured by a typical observer can be directly applied to the 
models discussed here. Changing from the coordinates (r, 9, (p) to a new set (r', i^', 9') such that the point r = rp, 6* = 
(with ro ^ 1) is the new origin of coordinates, one finds that the perturbation of the field (j) around r' — Q can be 
described as 5(j) = (j)o{t){'-f /2)\n. f , where / = (coshr' -|- sinh r'cos9') and tpo is the value of the field at the point 
r' — 0. The corresponding gauge invariant potential at horizon crossing is 

3 HtS^P 



The effect on the microwave background temperature fluctuations can be computed by integrating the Sachs- Wolfe 
effect along the line of sight ||2^ . The dominant effect is in the quadrupolc |jl^ , and it is of order 



5 mPp 2 



(51) 



Subcurvature fluctuations cannot by themselves give rise to inflating islands since their size is smaller than the curvature 
scale. 
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ST 
~T 



^10"'^?(l-^^o). (52) 



This is just a very rough order of magnitude estimate, which works well for JIq ^ 0.3. A more detailed study of the 
power spectrum of temperature anisotropics will be presented elsewhere psf . 

Note that if the universe is sufficiently flat, the factor (1 — 51o) rnay completely erase the effect. Otherwise, for a 
universe with appreciable curvature, we obtain a constraint on 7, 

7<7r^^lO-. (53) 



(1 - f^o) 



where we have used 5T/T < 10^*^, from the requirement that this effect does not dominate the temperature 
anisotropics on large scales, as seen by COBE [Q. Since Ht is larger than m^, it is clear that 7 has to be very small 
in order to avoid large temperature fluctuations. 

In the supernatural model, 7 ^ i?gm^/2, this constraint implies 

Ro<2x 10^2(1 - noy^^^H^'^ 

which is not difficult to accomodate. Note that the size of the bubble Rq is necessarily less than the Hubble radius 
in the false vacuum and hence can easily be much less than the Hubble radius in the true vacuum. However, in the 
decoupled model discussed in |jl^ the constraint (^sh forces Ht to be quite large, and this causes a problem of large 
quantum fluctuations in the supercurvature modes Q . 

In the class of models (|39|), the size of 7 is determined by Eq. (|4|). Clearly, the effect can be made small by 
choosing parameters such that mp and the size of the bubble are sufficiently small. This is not always straightforward 
to implement. For instance, the "hybrid" open inflation model considered in Ref. turns out to be quasi-open and 
suffers from too large semiclassical anisotropics. However, it is possible to write an open hybrid model, with a massless 
inflaton in the false vacuum, that satisfies the constraints |2q]. 



B. Supercurvature anisotropies 

In the previous subsection we have considered the case where the I = mode was "oversized" , meaning that it 
took an amplitude much larger than its expected r.m.s. Because of this, an observer far from the center of the 
inflating region would see the anisotropy (pVi). In this section we shall estimate the anisotropies caused by the I > 
supercurvature modes. Here we are not thinking that these higher modes are "oversized" ; they simply take random 
values of the order of their r.m.s. For simplicity, we shall consider an observer located at r = 0, but the effect should 
not be much different for an observer located elsewhere. 

The size of CMB anisotropies caused by the I > supercurvature modes has been estimated in . For the class 
of models (|3^), where the supercurvature mode is normalized as in (^4|), the quadrupole CMB anisotropies are of 
order 



ST 
~T 



sup 



1=2 ^ °'Ht T 



sub 



1^2 



(54) 



Here {ST/Ty^^ are the temperature anisotropies caused by the subcurvature modes (with > 0). The supercurvatre 
effect decreases very fast with multipole number, basically as (1 — J7o)'^^. If the fluctuations we observe in the CMB 
are due to inflation, then we need (5T/T)*"'' ~ 10^^, and from ( |5^ ) we have that Hp cannot be too much larger than 
Ht: unless the universe is almost flat. 

For the supernatural model, the supercurvature mode ( |28| ) has a normalization 2/{HpRo) times larger than its 
counterpart (M) [we are ignoring the mild enhancement due to the factor (w/df,)]. Hence, the analog of (p4) is 



ST_ 



sup 



1^2 



(i_f7„)^^ 

^ 'HtRo T 



sub 



1^2 



(55) 



Therefore we need Rq ^ H^. Since Rq has to be necessarily smaller than Hp^, we have a two- fold restriction. On 
one hand, Rq ^ Hp, and on the other Ht ~ Hp. Thus, it seems fair to say that the model is not as natural as 
it was thought to be jsj: the difference in energy density between the true and the false vacuum cannot span many 
orders of magnitude. The reason is the following: In spite of the fact that the fleld is massive in the false vacuum, a 
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supercurvature mode exists. Its normalization is not proportional to Hp as in the usual case (|39[), but to , which 
is even larger. The effect can be thought as the excitation of the pseudo-Goldstone modes due to the acceleration of 
the domain wall "boundary" . The model may still be viable in a certain range of parameters. Determining this range 
requires detailed analysis, which is left for future research 



VII. CONCLUSIONS 



Open inflation is an appealing way of reconciling an infinite open universe with the inflationary paradigm. In 
this scenario, a symmetric bubble nucleates in de Sitter space, and its interior undergoes a second stage of slow-roll 
inflation to almost flatness. Single-fleld models of open inflation can in principle be constructed, but it does not seem 
possible to do so without a certain amount of fine tuning 1^ . The basic problem is that there is a hierarchy between 
the large mass needed for successful tunneling and the small mass required for successful slow-roll. For that reason, 
it seems natural to consider two-field models of open inflation where one field does the tunneling and the other 
drives slow-roll inflation inside the bubble. 

In this paper we show that a large class of two-field models of open inflation do not lead to infinite open universes, 
as it was previous thought, but to an ensemble of inflating islands of finite size. The reason is that the quantum 
tunneling does not occur simultaneously along both field directions, and the equal-time hypersurfaces in the open 
universe are not synchronized with equal-density or fixed- field hypersurfaces. Technically, one finds that there are no 
0(4) invariant instantons for the two-field system which would describe the formation of a bubble with "large" values 
of the slow-roll field in its interior. Large values of the inflaton field, needed for the second period of inflation inside 
the bubble, only arise as localized fluctuations. The interior of each nucleated bubble will contain an infinite number 
of such inflating regions, giving rise to a rather unexpected form of the large scale structure of the universe in these 
models. 

The picture is the following. Right after the bubble has nucleated there will be, on the t = const, hypersurfaces 
inside the bubble, a certain density of occasional large fluctuations of the slow-roll held that lead to inflating islands. 
This fluctuations are caused by modes whose wavelength is larger than the curvature scale. Denoting by 7 the 
eigenvalue of the Laplacian on the unit Hyperboloid (with 7 < 1), the co-moving size of the inflating islands is given 



by d ^ 7-1 [The parameter 7 can be determined in terms of the parameters of the model, see Eqs. ( pO[ ) and (48), and 
it is important in deriving observational constraints]. Each one of the inflating islands will be a quasi-open universe. 
Since the volume of the hyperboloid is infinite, inflating islands with all possible values of the held at their center will 
be realized inside of a single bubble. We may happen to live in one of those patches where the universe appears to be 
open. The fact that the inflating regions are finite gives rise to classical anisotropics like those discussed in Ref. . 

In particular, we have studied the supernatural model introduced by Linde and Mezhlumian ||^. We have shown 
that in spite of the large mass of the inflaton fleld in the false vaccuum, there is a supercurvature mode. Its amplitude 
is proportional to R^^ ^ rather than the usual Hp. Here i?o is the radius of the bubble at the time of nucleation and 
Hp is the Hubble rate in the false vacuum. Since R^^ > Hp, this effect is quite important. In order to make the 
model compatible with observations, it is required that the energy density in the false vacuum should not be much 
larger than in the true vacuum. This means that Hp /Hp cannot span many orders of magnitude, as it was previously 
believed The supercurvature mode can be understood as the pseudo-Goldstone mode associated with the choice 
of a tunneling direction in field space. Combining the supercurvature anisotropies with the classical ones we find that 
the range of f2o will also be restricted. Detailed analysis is required in order to determine the range of parameters in 
which the model may still be viable [ p8[ . 

For the more general class of models (|9|), the size of the infiating islands can be chosen to be comfortably large 
by an appropriate choice of parameters. In this way, the classical anisotropy will be unobservably small. By order 
of magnitude, the constraint is given by Eq. (p3|), where 7 is given in Eq. (^). The constraint will be satisfied if 
the mass of the slow- roll field is sufficiently small in the false vacuum and Rq is much smaller than Hp^ . In a future 
publication p8[ we will give more precise constraints from the observed power spectrum of temperature anisotropies 
of the CMB. 

Finally, there are some two-field models of open inflation, such as the one introduced by Green and Liddle |^ in 
the context of induced gravity, which need not be affected in principle by the classical anisotropies mentioned above. 
In these models, the value of is not variable; it is determined in terms of the parameters in the potential. It would 
be interesting to check whether 0(4) symmetric instantons do indeed exist in this model. 
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APPENDIX A: 

We will show in this Appendix that the two-point function on a i = const surface for the state ip~' dies off as e"''''^/^, 
where d is the comoving distance between the points. We will compute the two-point function ouside the lightcone, 
and then continue it to the inside. 

To compute the two-point function for Lp2, hereafter Gi^j (x, x'), we will use the fact that the modes ypim are 
properly normalized as Klein-Gordon modes of mass + \ living in the r = const. (2-1-1) de Sitter hypersurfaces of 
the outside lightcone metric. Thus we can define the fields 



and the de Sitter invariant vaccuum anihilated by "^^a^j^. Notice that (Elds in ^ corresponds to the line element 
of a closed coordenization of a (2+1) dc Sitter space, for which the two-point functions can be found in Ref. ||3l| ]. 
We can now write the two-point function G^p^ in terms of the two-point functions for 3^p, 

pirn 

= E Fp[T)l]J?)dsmp{x')yp[x''mds. (A2) 
p 

where the sum over p has to be understood like a sum over the discrete eigenvalues of the Schrodinger equation ( p^ ) 
and like an integration over its continuum spectrum > 0. On a given r — const hypersurface, the p dependence 
of {x,x') will be given by Gj,(x% x'*), weighted for each p by a^j^l |-Fp| p. The two-point function Gp can be found 
in ||] 

r f'\ 1 r(l-^p)r(l + ^p) ^^^ , . . 3 1 + Z\ 

where F is the hypergeometric function and Z is the scalar product of the position vectors at points x^ and x'* in the 
embedding (3-1-1) Mikowski space, 

Z (x* , x"- ) = ^'^ (x* )^^ (x'* ) = cos J cosh p cosh p' ~ sinh p sinh p' . ( A4) 

Here 7 is the angle on the 2-sphere between the two points. We recall that for the lowest discrete eigenmode to first 
order in the shift 7, = 1 — 7/2, so we will denote by G'^ the two-point function for this eigenmode. 

Now we have to analytically continue ( |A3| ) to the inside of the lightcone by means of (|^). This amounts only to 
analytically continuing the scalar product Z, 

Z(x\x'*) —^ — C0S7 sinh r sinh r' + coshr coshr'. (A5) 

Taking r' = and r — so Z — coshd, and using eq. (9.131.1) in Ref. |32j, we find that inside the light-cone the 
two-point function between points separated a comoving distance d can be written as 

1 r(2-7/2)r(7/2) ^ l-coshd V^/^^^7 7-I 3 l + coshd\ 
^W-j^ ^ ) ^(^2'^'2'coshd-lj- 



As d — > 00, the hypergeometric function in (A6) tends to a constant, and the assymptotic behaviour of G'^ is given by 

1 r(7/2)r(2-i) /-iV"'\-,i/2 



which dies off exponentially with d. 

Here, we have only computed the first term in the sum (A2). The terms with > decay as e"*"^^, and hence 
they are subdominant at large distance. 
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APPENDIX B: 



To compute {/S.^IY for large r, we will need the asymptotic expressions for the hyperbolic harmonics 3^A,im for 
r ^ 7^^. The Legcndre functions are given by (see e.g. Ref. [ ^ ) 



p-/;,(cosh.) = ^-^j^, 

From equation the supercurvature mode 3^a,oo is given by 



where A = 1 — 7/2. For large r, we have 



1 



r(i-A)r(i + A) 



1/2 



sinh Ar 
A sinh r ' 



g-r7/2 

:V(l-7/2),00^ ^^[1 + 0(7)]- 



7rV2^ 

For Z > 0, we can express Pi, ' in terms of P^^^l^ using the recursion formula 

— /— 1/2 

The Legendre function Pu acquires then the form 

"sinh'r ^ ' 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



k=-l 



where Cfe(i^) are some functions depending on u. In fact, for large r, we do not need to compute all C'k{i'). We have 

t o t ake into account that for a supercurvature mode, P,^ (coshr) behaves for large r as e'"''' (as can be seen from 
(Bl)). Thus, for the supercurvature mode v = (1 — 7)/2, the term A: = Z in (B5) grows exponentially faster than the 
rest of terms in the sum , so the main contribution for large r will be given by this term. The coefficienty Ci{v) can 
be easily read from (B4): 



CM 



T{v + l/2) 



For large r, using 

3^(1-7/2),;^ 



and (23), we obtain 

r(Z + 7/2) 



1/2 



r(l-7/2)e-W2^™(f]). 



_^(l-7/2 + 0(^-7/2)r(Z-7/2). 
Finally, using \\mx-*oo r(x + a)/T{x) — x"", we can write 3^(i_-y/2),im for large Z, to order 7, as 

;7/2-l 



3^(1-7/2), in 



-'•^/2l^™(17)[l + 0(7)]. 



Using the results derived above, we can compute the amplitude of the I — mode near the origin, 

2 



1 



-] [1 + 0(7)], 



and for large r. 



7 \2 



(¥'7=0 



1 / 



(^7=0)" 



r-O 



(B6) 



(B7) 



(B8) 



(B9) 



(BIO) 
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As we can see, the amplitude of the mode I = decays exponentially for r ^ 7"^. Taking into account that we 
have choosen an 0(3, 1) symmetric vaccuum, this decrease in amplitude for large r must be compensated by the joint 
contribution of the I > modes, smeared over a suitable length scale, in such a way that the r.m.s fluctuations of the 
field are independent of r. Let us check it. We need to compute 



1=1 m=-l 1=1 

where is a certain cutoff, which has to grow as we move away from the origin to include more and more modes in 
the sum. If we smear the field over a fixed comoving lenght ^, realizing that the wavelength of the ^-th multipole is 
proportional to (sinhr)/Z, we can take h = sinhr/^. Finally, we obtain 

1 / CTo 



(A^lr-^l^^ j(20-"- (Co)1..o(l-7ln2e). (B12) 

As we can see, as long as | ln2^| ^ 7^^, the added contribution of the relevant modes is the same as the one given by 
the / = mode near the origin. 



APPENDIX C: 

In the thin wall approximation, neglecting gravitational backreaction, the background geometry is found ||25|| to be 
described by two de Sitter pieces with different Hubble constant glued together at some rjw ■ The scale factor is given 
by 

ttEiv) = apivWiv - Vw) + arivWivw - v), (CI) 
where ap and ut are the scale factors in the false and in the true vaccuum, 

1 



Hp cosh 77 

AT (77) 



Ht cosh(?7 — S) 
Continuity of a^; at the wall implies 



ainw) = 77 \ = 77 ^7 K = ^0, (C2) 

Hpcoshriw Ht cosh(riw — 0) 



where Rq is the radius of the wall, and S is given by 



e 



(1 + ^l-H^Rl)il ^l-HlRl) 



HpHtRq 

To complete the description, we need to know the value of Rq. It can be found in Ref. 



(C3) 



i?o - ""^^ ^ (C4) 

^{Hl -H'^ + (k5i/2)2)2 + Km'iSf 

where k = SnG and 5*1 is the wall tension. 

We want to find the lowest eigenvalue of the Schrodinger equation ( p^ in the background given above. The effective 
potential is given in this case by 

U = a|[m2 + ga^ _ 2(ff| 9{rj - rjw) + 0{r^w - v))] + (^F - nT)S{v ~ Vw), (C5) 

where Hp — a'p/ap and similarly for Ht- 
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We will take a perturbative approach. We will divide the effective potential lA in an unperturbed one, IAq, and in a 
small perturbation part, }JAi: 

Uq = -iapHl (C6) 
}Ux = al (m2 + gal - 2H'^e{riw - v)) + (^f - nT)S{r] - Vw) 

+2Hlal0{r,w -V)- (C7) 

The unperturbed Uq corresponds to the effective potential of a massless scalar field in de Sitter space, which has as a 
ground state a supercurvature mode with energy Pq = —1 and wavefunction |^ 

F_i = ^Mv)- (C8) 

To first order in perturbation theory, the shift of the energy Pq = — 1 is given by 

7 = - 1) = + ^{m^^ - mil (C9) 

where ttif is the efective mass of the slow- roll field in the false vaccuum, and rriT the effective mass in the true 
vaccuum. In this case, m?p = m? and m3p = m? + gv^ . 
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